Nuclear pygmy modes and the dynamics of the nuclear skin 
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The information on pygmy resonances reveals new aspects on the isospin dynamics of the nucleus 
with important astrophysical consequences. In this connection, the precise knowledge of nuclear 
response functions plays a key role in the determination of photonuclear reactions cross sections 
which are of importance for the synthesis of heavy neutron-rich elements. For that purpose, a 
theoretical method based on density functional theory and multi-phonon approach is applied for 
investigations of nuclear excitations with different multipolarities and energies in stable and exotic 
nuclei. The possible relation of low-energy modes to the properties of neutron or proton skins is 
systematically investigated for isotonic and isotopic chains. Our studies of dipole and quadrupole 
response functions and the corresponding transition densities indicate new pygmy dipole and pygmy 
quadrupole resonances, describing oscillations of the nuclear skin. Also, the presence of skins is found 
to affect the magnetic response of nuclei. 



I. INTRODUCTION 

The studies of nuclear electromagnetic response along 
isotonic and isotopic chains of nuclei led to surprising new 
observations of new modes of excitations different from 
known surface vibrations and giant resonances in nuclei 
close to the isospin N=Z symmetry. A clear example is 
the existence of an enhanced dipole strength located be- 
low or closely above the neutron particle emission thresh- 
old which was found to be a common feature of stable and 
unstable nuclei with neutron excess HHl. As far as this 
bunching of dipole states with very similar spectroscopic 
features resemble a resonance, the new mode was named 
pygmy dipole resonance (PDR). Even though, the na- 
ture of the PDR is a subject of outstanding discussions, 
the most common theoretical interpretation associated it 
with oscillations of a weakly bound neutrons (for N/Z>1) 
or protons (for N/Zsal) with respect to an isospin sym- 
metric core • A question, coming up immediately in 
this connection, is to what extent the presence of a neu- 
tron or proton skin will affect excitations of other mul- 
tipolarities and parities and vice versa. Promising can- 
didates are low-energy 1 + , 2 + , 3 _ states ect., especially 
those in excess of the spectral distributions known from 
stable nuclei. Recently, quadrupole response functions 
were investigated, theoretically in neutron-rich Sn nuclei. 
A close connection of low-energy 2 + excitations and nu- 
clear skins was found. These quadrupole states were re- 
lated to pygmy quadrupole resonance (PQR) [l(| ■ In this 
connection, QRPA and QPM calculations of dipole and 
quadrupole strength distributions for pygmy resonances 
in nuclei with Z=12, 13,50 and N=50,82 are discussed. 
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H sp + Hp a i r is the mean-field part defining 
the single particle properties including pairing interac- 
tions for protons and neutrons. Different from the stan- 
dard QPM scheme this part is obtained self-consistently 
by a fully microscopic Hartree-Fock-Bogoliubov (HFB) 
approach (l2l |. H^f, H^ M and are residual inter- 
actions, taken as a sum of isoscalar and isovector sep- 
arable multipole and spin-multipole interactions in the 
particle-hole (ph) and multipole pairing interaction in 
the particle-particle (pp) channels. The model param- 
eters are fixed either empirically 13] or by reference to 
QRPA calculation performed within the density matrix 
expansion (DME) of G-matrix interaction discussed in 
Ref. [lj. 



A. The nuclear ground state 

The reliable description of ground state properties is 
of genuine importance for extrapolations of QRPA and 
QPM calculations into unknown mass re gion s. Taking 
advantage of the Kohn-Sham theorem [l4l,|l5j of density 
functional theory (DFT) the total binding energy B(A) 
of the nucleus could be expressed as an integral over an 
energy density functional with (quantal) kinetic (r) and 
self-energy parts, respectively, 



B(A) = J d 3 
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II. THE THEORETICAL MODEL 



J pair 



The model Hamiltonian resembles in structure the 
standard quasiparticle-phonon model (QPM) [11] but in and pairing contributions are indicated by E pa i r . The 
detail differs in the physical content in important aspects second relation is obtained from HFB theory with occu- 
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FIG. 1. BCS ground state densities of Sn isotopes obtained by 
the phenomenological DFT approach and used in the QPM 
calculations. 



pancies u| and potential energies < 



U >j of the occu- 



pied levels j, see e.g. [12|. Above, U(p) is the proper 
self-energy, i.e. not including the rearrangement contri- 
butions from the intrinsic density dependence of nuclear 
interactions [12j. Hence, U(p) has to be distinguished 



from the effective self-energy obtained by variation 

1 dppU(p) 

2 dp 
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and appearing in the single particle Schroedinger equa- 
tion. In order to keep the QPM calculations feasible we 
choose S = Uws to be of Wood-Saxon (WS) shape with 
adjustable parameters. By inversion and observing that 
the densities and potentials in a finite nucleus are nat- 
urally given parametrically as functions of the radius r, 
we find 



p{r)U(r) = -2 



ds^-U ws (s) . (2.4) 
as 



Evaluating these relations with the microscopic pro- 
ton and neutron densities obtained by solving the 
Schroedinger equation with Uws the potential U(p) is 
the self-consistently derived reduced self-energy entering 
e.g. into the binding energy. 

In practice, for a given nucleus of mass A the depth 
of the central and spin-orbit potentials, radius and dif- 
fusivity parameters of Uws are adjusted separately for 
protons and neutrons to the corresponding single par- 
ticle separation energies, the total binding energy [16j . 
the charge radii and (relative) differences of proton and 
neutron root-mean-square (RMS) radii, 
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FIG. 2. BCS ground state densities in N=50 obtained by 
the phenomenological DFT approach and used in the QPM 
calculations. 



from our previous HFB calculations [12]. In Ref. [8| 
theoretically obtained RMS radii are compared to those 
determined from charge exchange reactions for a number 
of Sn isotopes. The approach sketched above leads to 
very satisfactory results on binding energies and proton- 
neutron RMS-differences as shown in 0, Q • A smooth 
dependence of the parameters on A is found which sup- 
ports the reliability of the method. 

Calculations of ground state neutron and proton den- 
sities for Z=50 and N=50 are shown in Fig. Q] and Fig. 
[U Results from our calculations in N=82 isotones could 
be found in Ref. [||. Of special importance are the nu- 
clear surface regions, where the formation of a skin takes 
place. A common observation found in the investigated 
isotopic and isotonic chains of nuclei is that the thickness 
of the neutron skin is related to the neutron-to-proton ra- 
tio N/Z. Thus, in Sn isotopes with A>106 the neutron 
distributions begin to extend beyond the proton density 
and the effect continues to increase with the neutron ex- 
cess, up to 132 Sn. Thus, these nuclei have a neutron skin. 
The situation reverses in 100 ~ 102 5n, where a tiny proton 
skin appears. For the case of N=50 isotones the neutron 
skin thickness decreases from 86 Kr toward 92 Mo which 
has the largest proton number Z. 



B. The nuclear excited states 

The nuclear excitations are expressed in terms 
of quasiparticle-random-phase- approximation (QRPA) 
phonons 
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where j = (nljmr) is a single-particle proton or neu- 
tron state; and are time-forward and time- 
backward operators, coupling two-quasiparticle (2QP) 
creation or annihilation operators to a total angular mo- 
mentum A with projection fx. The excitation energies 
of the phonons and the time-forward and time-backward 
amplitudes and <fj*j 2 in Eq. (I2.6[) are determined 
by solving QRPA equations [llj ]. 

Furthermore, the QPM provides a microscopic ap- 
proach to multiconfiguration mixing [11]. For spherical 
even-even nuclei the model Hamiltonian is diagonalized 
on an orthonormal set of wave functions constructed from 
one-, two- and three-phonon configurations [l7| . 
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where R, P and T are unknown amplitudes, and v labels 
the number of the excited states. 

The electromagnetic transition matrix elements are 
calculated for transition operators including the interac- 
tion of quasiparticles and phonons [l8j where exact com- 
mutation relations are implemented which is a necessary 
condition in order to satisfy the Pauli principle. 

The QPM allows for sufficiently large configuration 
spaces such that a unified description of low-energy 
single- and multi-particle states is achieved also for the 
GDR energy regions. Such a unified treatment is exactly 
what is required in order to separate the multi-phonon 
and the genuine PDR 1~ strengths in a meaningful way. 



III. DISCUSSION 

From QRPA calculations in 110_132 Sn, a sequence of 
low- lying one-phonon 1 — states of almost pure neutron 
structure at excitation energies E* = 6 — 7.5 MeV is 
observed. These states correspond to excitations of least 
bound neutrons from the 3s-, 2p- and 2d-subshells with 
only a minor proton contribution of less than 1%. Similar 
results are found also in N=50, 82 isotones [!, H|. These 
low-energy dipole states have been related to neutron 
PDR. 

In the lightest 100 ~ 104 Sn isotopes, the lowest dipole ex- 
citations at E*= 8.1-8.3 MeV are dominated by proton 
excitations involving quasibound 2p 3 / 2 and lc/9/2 proton 
states [8j . Electromagnetic breaking of isospin symmetry 
is the main reason for the persisting of low-energy dipole 



strength close to 100 Sn. QRPA calculations of dipole re- 
sponse functions in several tin isotopes are presented in 
Fig. [31 The observed sequence of low-energy neutron or 
proton dipole states have been related to neutron or pro- 
ton PDR, respectively Q. For a more detailed insight 
into the characteristic features of the dipole excitations 
we consider the evolution of the proton and neutron tran- 
sition densities in various energy regions. In Fig. U (left 
panel) we display dipole QRPA transition densities for 
several tin isotopes for three different regions of excita- 
tion energies: the low-energy PDR region below the neu- 
tron emission threshold, the transitional region up to the 
GDR and in the GDR region and beyond. The transition 
densities displayed in Fig. 2] (left panel) were obtained 
by summing over the transition densities of the individ- 
ual one-phonon 1~ states located in the energy intervals 
denoted at the top of each column of the figures. A com- 
mon features of the presented 112 > 122 . 1 32g n cases j s fag^ 
up to E* = 8 MeV the protons and neutrons oscillate 
in phase in the nuclear interior, while at the surface only 
neutron transitions contribute. This pattern is generic to 
the lowest dipole states making it meaningful to distin- 
guish these excitations from the well known GDR states. 
Hence, we are allowed to identify the PDR as a distinct 
and unique excitation different from the giant dipole res- 
onance (GDR). Similar results for the QRPA neutron 
and proton dipole transitions densities are obtained also 
in N=50 isotopes which are presented in Fig. [4] (right 
panel). The exited QRPA 1~ states located below exci- 
tation energy E* <9 MeV are related to PDR. 

The transition of the neutron PDR to a proton PDR 
manifests itself via the proton and neutron transition 
densities as well. Thus, in 100 Sn the proton oscillations 
dominate at the nuclear surface due to the formation of a 
proton skin (see in Ref. Q). In comparison, the presence 
of a neutron skin explains the neutron oscillations at the 
surface of tin nuclei with A<106 shown in Fig. [T] (left 
panel) . 

In Fig. [5] the total PDR strength calculated in Z=50 
(left panel) and N=50 (right panel) nuclei is presented in 
comparison with the neutron skin thickness Sr, eq. (|2.5[) 
of these nuclides. It is shown that the two quantities 
are closely related. Here, the B(E1) sums are taken over 
QRPA dipole states associated with PDR according the 
analysis of transition densities and the structure of the 
state's vectors, respectively (see in Ref. [§[ for more de- 
tails). The close relationship of the PDR strength and 
the skin thickness is also confirmed in our investigations 
of N=82 nuclei @,[|. 

The comparison of QRPA and two- and three-phonon 
QPM calculations of dipole strengths indicates that for 
the PDR region the coupling of QRPA PDR and GDR 
phonons and multiphonon states is very important. The 
result is a shift of El strength toward lower energies 
which is well visible in Fig. [6] (left panel) . Consequently 
in the energy interval < E < 8 MeV the QPM calcu- 
lations give about twice as much total El-strength than 
the QRPA calculations, as shown in Fig. [B] (right panel). 
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FIG. 3. QRPA calculations of electric dipole strength distributions in Sn isotopes. 




FIG. 4. Dipole one-phonon proton and neutron transition densities in (left panel) Z=50 isotopes and (right panel) N=50 
isotones. 



Differently, with the increase of the excitation energy 
toward the GDR the QRPA and QPM summed B(E1) 
transition strengths are of comparable size which sug- 
gests the prevalence of p — h excitations of GDR type. 
Thus, our QPM calculations of 112 Sn (N/Z=1.24) and 
120 Sn (N/Z=1.4) confirm experimental and QRPA ob- 
servations of increasing PDR strength with the increase 
of N/Z ratios 0. The comparison of 112 Sn and 120 Sn nu- 
clei shown in Fig. [B] (right panel) shows that the relative 
difference between the El strengths calculated for certain 
energy regions, both for the QRPA and QPM cases, re- 
duces with the increase of the excitation energies toward 
GDR. This behavior is expected because of the weak de- 
pendence of GDR on N/Z ratios in tin isotopes. 

Furthermore, in theoretical investigations of 2 + exci- 
tations in Sn isotopic chain, we find a strength clustering 
of quadrupole states, at low-energies, with predominantly 
neutron structure. At the same time, the proton contri- 



bution to state vectors and B(E2) strengths, located in 
the energy range E* « 2-4 MeV, increases toward 104 Sn 
and brings to more intensive proton quadrupole excita- 
tions in 104 Sn there. Consequently, the 104 Sn nucleus 
appears to be an opposing case where a change from a 
neutron to a proton skin occurs. Theoretical results of 
B(E2) strength distributions in 104 420,i34g n j so ^ p es are 

presented in Fig. [7J A sizable increase of B(E2) strength 
at E* ps 2-4 MeV is observed for the heaviest tin isotopes 
- 130 Sn and 134 Sn - studied in 10]. The clustering of 
quadrupole states at low-energies shows a pattern simi- 
lar to the PDR phenomenon. Therefore, we may consider 
the spectral distribution of a Pygmy Quadrupole Res- 
onance (PQR). Correspondingly, the connection of the 
PQR with neutron or proton skin oscillations is demon- 
strated in analysis of transition densities (see in Ref . [13] ) ■ 
Similarly to the PDR, a transition from a neutron PQR 
to a proton PQR in 1 Sn is observed for the mass region 
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FIG. 5. (left panel) QRPA results for the total PDR strengths in the 100_134 Sn isotopes are displayed for comparison together 
with the nuclear skin thickness 8r, eq. (|2.5|l . Experimental data on the total PDR strengths in 112 . 116 > 120 > 124 Sn [19,20]) are 
also shown. In the lower panel, the skin thickness derived from charge exchange reactions by Krasznahorkay et al. published 
in Refs. [21] are indicated, (right panel) QRPA results for the total PDR strengths in N=50 isotones. 




FIG. 6. (left panel) Calculations of electric dipole strength up to E*= 9 MeV in 120 Sn: a/QRPA, b/ two-phonon QPM and c/ 
three-phonon QPM; (right panel) Comparison of QRPA and three-phonon QPM B(E1) strengths summed up to E*= 8 MeV 
(PDR region), 8.5 MeV (GDR region) and 9 MeV (GDR region), respectively. 



where the neutron skin reverses into a proton skin [loj 
and in Fig. [5] (left panel). 

From the detailed analysis of the state vectors struc- 
ture of the 2 + excitations in Sn nuclei at excitation ener- 
gies E* = 2 — AMeV we find that it is dominated by neu- 



tron 2QP excitations from the valence shells. The most 
important proton contribution in all isotopes is due to 
the [139/22^5/2]^ 2QP component, which, however, never 
exceeds 5% [Kj. This dependence is illustrated for the 
[^t]QRPA states in Fig. |8] (left panel). The figure re- 
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FIG. 7. QRPA calculations of electric quadrupole strength distributions in Sn isotopes. 



99,9 
99 
95 

/-s 80 
^ 60 
w 40 

3 20 
5 
1 

0,1 



(2d 



■ ■ " ■ 

(2d„x3s,J 

1 s,! " 2 - (2d x3s ) 



neutrons ) 
protons 



<2d r* i^^.. 

>'* (21.,,/. ■ 



, (2d SI! x3s 1J J - 
■ 

< 1 9») ! . 

■ 



< lh „ 



< 2d »2 x3s „2>, (1g ra x2d I(! ), 



'9/2 5/2 'v. 
• • • • 



100 105 110 115 120 125 130 135 

A 




1E-3 



100 105 110 115 120 125 130 135 

A 



FIG. 8. (left) QRPA structure of the [I^Iqrpa states calculated in Sn isotopes. The main neutron and proton 2QP components 
with a state vector probability Wj 1 j 2 (Ai)=|^ 1 ! 3 - 2 | 2 — Iv'jjjJ 2 larger than 0.1% are presented; (right) QRPA results for B(E2)f 
transition probabilities summed over the [2%]qrpa and [2^}qrpa states in Sn isotopes. 



fleets the change of the proton binding ener gy e b of the 
(?9 j2 level when approaching the N — Z limit [10( . In the 
neutron sector, the contributions follow closely the evolu- 
tion of the shell structure as seen from Fig. [S](left panel). 
In most cases the ft^QRPA state vectors are dominated 
by re-scattering contributions related to re-orientation of 
the s.p. angular momenta [To| . 

QRPA calculations of summed B(E2) strengths in 
i04-=-i34g n isotopes are presented in Fig. [5] (right panel). 
A sizable increase of B(E2) values summed at E* ss 2-4 
MeV is observed with increasing of N/Z ratio [13 as- 
sociated with the mass dependence of the neutron 2QP 
transition matrix elements and decreasing neutron bind- 
ing energy e&. For N/Z < 1.4 the summed B(E2) values 
follow the increasing contribution of the protons, which 
are coupled directly by their physical charge to the elec- 
tromagnetic field. 



An effect of increased low-energy E2 strengths with 
increasing neutron number is observed as well in our 
QRPA calculations of lighter mass nuclei 33 A1/ 35 A1 and 
32 Mg/ 34 Mg which are presented in Fig. |H] (left and right 
panels), respectively. That mass region is of particular 
interest because of the well known 'island of inversion' at 
the N=20 shell closure 

Correspondingly, the connection of the low-energy 2 + 
states with a PQR is demonstrated in the analysis of 
transition densities [l(J ■ Results from QRPA proton and 
neutron quadrupole transition densities of the PQR mode 
in 124 Sn are presented in Fig. [TU](left panel). Strong neu- 
tron oscillations at the nuclear surface play a dominant 
role in the PQR energy range which is in agreement with 
Ref. 0. Similarly to the PDR, a transition from a 
neutron PQR to a proton PQR in 104 Sn is found for the 
mass region where the neutron skin reverses into a proton 
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FIG. 9. QRPA calculations of B(E2) spectral distributions normalized to EWSR: S n = KE?B(E2) in (left panel) 33 Al and 
35 Al and (right panel) 32 Mg and 34 Mg. 
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FIG. 10. (left panel)QRPA calculations of quadrupole transition densities related to the main proton and neutron components 
in the state vectors of [2f\QRPA states with excitation energies E*= 2-5 MeV in 124 Sn; (right panel) Multiphonon calculations 
of B(E2) transitions probabilities of low-energy 2 + excitations in 124 Sn. 



skin [T(j. Furthermore, QPM calculations of B(E2) and 
B(M1) transition rates of low-lying 2 + states in 124 Sn 
indicate clearly differences from known collective states 
and scissors modes. The results on B(E2;g.s — > 2+) 
transition probabilities are presented in Fig. [10] (right 
panel)). 

An important question to clarify is the fine structure 
of the observed low-energy dipole strengths. For this 
purpose, QPM calculations of low-energy El and spin- 
flip Ml excitations are made, in a configuration space 
including up to three-phonon components, in 138 Ba nu- 
cleus. The results are compared to experimental data 
obtained from polarized photon beams |5j. The experi- 
mental observations and their theoretical interpretation 



show unambiguously the predominantly electric charac- 
ter of the observed low-energy dipole strength (see Fig. 
HU). 

In conclusion, in systematic investigations of low- 
energy excitations of different multipolarity in Z=12, 13 
, 50 isotopic and N=50, 82 isotonic chains a specific sig- 
nals of new modes of excitations related to PDR and 
PQR are observed. As common feature of the structure 
of the PDR excited states in neutron-rich nuclei is dom- 
inated by neutron components while with the increase 
of the proton number, toward N=Z it transforms into a 
proton one. Consequently, the PDR transition strength 
is found directly related to the thickness of a neutron or 
proton skin. Its generic character is further confirmed 
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FIG. 11. Experimental data (left panel) and QPM calculations (right panel) for El (electric dipole) and Ml (magnetic dipole) 
photoabsorption cross sections in 138 Ba below 9.0 MeV. The cross sections are averaged over 0.2 MeV energy bins (see in Ref. 
5 as well). 



by related transition densities, showing that this mode is 
clearly distinguishable from the conventional GDR mode. 
In investigations of low-energy quadrupole states, PQR 
mode resembling the properties of the PDR is identified. 
An especially interesting aspect is the close relation of 
the PQR excitations to the shell structure of the nucleus. 
The most convincing evidence for this feature is the dis- 
appearance of the PQR component in the double magic 
132 Sn nucleus. Thus, the PQR states are containing im- 



portant information on the structure of the valence shells 
and their evolution with the nuclear mass number. 

In N=82 nuclei the fine structure of the low-energy 
dipole strength is explained as being mostly of electric 
character. Even though the magnetic contribution to 
the PDR is small in comparison with the electric one 
it should be considered of a great importance for the un- 
derstanding of the spin dynamics of the nucleus in the 
presence of skin. 
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